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Abstract 

We give the results of a study on the decay we performed in the Gran Sasso Laboratory (LNGS) by detecting the 
gamma rays from the radon progeny. The motivation was to monitor the stability of radioactivity measuring several 
times per year the half-life of a short lifetime (days) source instead of measuring over a long period the activity of 
a long lifetime (tens or hundreds of years) source. In particular, we give a possible reason of the large periodical 
fluctuations in the count rate of the gamma rays due to radon inside a closed canister which has been described in 
literature and which has been attributed to a possible influence of a component in the solar irradiation affecting the 
nuclear decay rates. We then provide the result of four half-life measurements we performed underground at LNGS in 
the period from May 2014 to January 2015 with radon diffused into olive oil. Briefly, we did not measure any change 
of the ^^^Rn half-life with a 810“^ precision. Linally, we provide the most precise value for the ^^^Rn half-life; 
3.82146(16)sto,(4)j3„, days. 
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1. Introduction 

A possible time dependence of the radioactive nuclei 
decay constant has been searched for since the begin¬ 
ning of the science of radioactivity. Lor instance, in the 
Ph.D. of M. Curie m one can already find the descrip¬ 
tion of the search for a difference in the radioactivity 
of uranium between midday and midnight. Recently, in 
particular since the year 2009 131, various experiments 
have reported evidence of a time modulation of the de¬ 
cay constant of several radioactive nuclei with period, 
in most cases, of one year (but also of about one month 
or one day) and amplitude at the per mil level. This an¬ 
nual modulation, with the maximum in Lebruary and the 
minimum in August, has been correlated to the change 
of the Sun-Earth distance between aphelion and perihe¬ 
lion; however different measurements exclude any mod¬ 
ulation as large as the reported ones (for recent reviews 
on results with and without modulation I[3]l4l). 

In 12 |5l the existence of new and unknown parti¬ 
cle interaction has been advocated to explain the yearly 
variation in the activity of radioactive sources, with the 


particle being emitted from the Sun. In Q the pos¬ 
sibility that the coupling to a long range scalar held, 
sourced by the Sun, might be the origin of the mod¬ 
ulation has been quantitatively discussed. The labora¬ 
tory constraints on the variation of aem and of the elec¬ 
tron to proton mass ratio on an annual timescale turned 
out to induce upper bounds to the relative variation of 
the decay constant nine orders of magnitude lower than 
the claimed per mil effect. Solar neutrinos have also 
been proposed as responsible for the modulation, with 
a cross section several orders of magnitude higher than 
expected. In particular, the possibility for anti-neutrinos 
affecting the decay of ^^Na has been very recently 
studied in a reactor experiment i). 

Since we believe that dedicated experiments are still 
needed to clarify the somehow contradictory situation, 
in this letter we describe the results of a different ap¬ 
proach we pursue to monitor the time dependence of 
radioactivity. 

Generally, a long lifetime source (tens or hundreds of 
years) is selected and its activity is measured for a pe¬ 
riod of at least one year. This way both the count rate 
and the dead time are almost constant. We also followed 
this approach in the study of ^^^Cs Q, and ^^^Th 
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Table I: The decay chain m. Gamma rays with relative prob¬ 
ability smaller than 3% are omitted. 


Isotope 

Decay 

type 

half-life 

Gamma 
energy [keV] 

Relative 

probability 


a 

3.8 d 



218po 

a 

3.1 m 



214pb 

p 

26.8 m 

242 

7% 




295 

18% 




352 

36% 

214Bi 

p 

19.9 m 

609 

45% 




768 

5% 




934 

3% 




1120 

15% 




1238 

6% 




1378 

4% 




1764 

15% 




2204 

5% 

214po 

a 

164^8 



210pb 

p 

223 y 

46.5 

4% 

210po 

a 

138 d 



206pb 






naiii. However, during such a long period of time sev¬ 
eral things may change, such as the laboratory pressure, 
temperature and humidity, the radioactive background, 
the electronic noise, and the performances of the elec¬ 
tronic chain. It is not easy to keep all these under control 
if looking for a variation at the per mil level, or even be¬ 
low. To mitigate the effect of all these possible instabil¬ 
ities we decided to measure time after time the half-life 
of a short lifetime (days) source and for this purpose we 
selected ^^^Rn, a source relatively easy to produce. A 
similar approach has already been followed in ||9l and 
IfTOll with the measurement of the '®^Au and ^''^Po half- 
life, respectively, with a precision of a few parts over 
lOl 

In the next section we describe the experiment we 
performed with radon in a confined air volume, before 
giving in section 3 the results of the measurements we 
performed underground at LNGS with radon diffused 
into olive oil. 

2. Radon in air 

^^^Rn is a noble radioactive gas coming from the 
decay of (half-life: 1600(7) y E] ) in the 

chain. The decay chain below ^^^Rn (half-life: 
3.8232(8) d IfT^ l is given in Table 1, together with 
the energy and the relative probability of the emitted 
gamma rays. In our experiment we measure the ^^^Rn 
half-life with gamma spectroscopy: this way we are sen¬ 
sitive to variations both in alpha and beta decay. 
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Figure 1: Normalized deviations from the exponential trend of the 
measured count rate as function of time for radon in air: without (up¬ 
per panel) and with (lower panel) the polystyrene particles inside the 
glass sphere containing air charged with radon. The shaded areas limit 
the region where fluctuations are within ±2cr. 

The set-up we used in the external Gran Sasso Lab¬ 
oratory during the preparatory phase of the experiment 
(to study the source production) is rather simple: a glass 
sphere (130 mm diameter) is connected through a pipe 
to a stainless steel cylinder containing a 0.3 kg rock rich 
in uranium. The radon from the radium decays in the 
rock fills by diffusion the glass sphere which, after 5- 
6 days, is isolated from the radon source by closing a 
Nupro (R) valve. 

Gamma rays from radon progeny are detected by a 
3”x3” Nal crystal placed at a few millimeter distance 
from the sphere surface (the geometrical centers of the 
sphere and of the crystal are aligned). The Nal is pow¬ 
ered and read out through an Ortec (R) digiBASE (TM) 
which also pre-amplify and digitize the signal. Both the 
detector and the glass sphere are enclosed inside a 5 cm 
thick lead shield. The detector count rate, i.e. the inte¬ 
gral of the spectrum above 6 keV threshold, is of about 
1000 counts per second (cps) at the beginning of the 
measurement, to be compared to a background, before 
filling with radon, of about 5 cps. 

The normalized residuals of the count rate, i.e. the 
difference between the measured rate and the expected 
one, divided by the expected rate as function of time 
is shown in Figure 1 (upper panel). Instead of having 
a statistical distribution around zero, there is clearly a 
24 hour period. We have the same behavior if, instead 
of considering the whole energy spectrum, we take into 
account only the peaks due to and ^'"^Bi decay. A 
similar effect has already been observed in the gamma 
radiation from radon progeny within air in confined con¬ 
ditions lEinKB] and it has been attributed to a pos- 
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Figure 2: Normalized difference of detector count rates (D1-D2) as 
function of time with the two detectors placed on opposite sides of 
the glass sphere filled with radon (top panel). Typical 2(t fluctuations 
amount to about 0.001. The time dependence of the temperature dif¬ 
ference (T2-T1) measured on the sphere surface is shown in the lower 
panel. The hours from 100 to 150 correspond to a week-end, when 
laboratory temperature is very stable. 

sible influence of a component in the solar irradiation 
affecting the nuclear decay rates. In our opinion, all this 
could simply be attributed to the displacement of the ra¬ 
dioactive nuclei inside the gas volume, with a variation 
in the detection efficiency. 

We tried to verify our hypothesis by completely Ail¬ 
ing the glass sphere with spherical polystyrene particles 
(diameter: 0.7-0.9 mm). This way the atoms cannot 
move freely within the entire volume but they are essen¬ 
tially confined to the small volume within the interstitial 
space among polystyrene particles. The background, 
measured during 20 days before Ailing with radon, had 
still a stable rate of about 5 cps. The glass sphere was 
then charged with radon for two weeks to give an initial 
count rate of about 1500 cps. Spectra taken with this 
configuration do not exhibit any modulation, as shown 
in Figure 1 (lower panel). As a consequence, we con¬ 
clude that the displacement of radioactive atoms is the 
reason of the effect. 

Searching for the trigger of this displacement, we 
modified the experimental set-up by adding a second 
3”x3” Nal detector (the two detectors were symmetri¬ 
cally placed on opposite side of the glass sphere). In ad¬ 
dition, we measured the temperature at different places 
on the sphere by using precision integrated-circuit tem¬ 
perature sensors with output voltage linearly propor¬ 
tional to the Centigrade temperature. In particular, we 
measured the temperature in the two points on the sur¬ 
face just facing the two detectors. 

We noticed that the temperature change is not the 


same everywhere on the surface of the glass sphere. In 
addition, we had a clear correlation (Figure 2) between 
the count rate difference of the two detectors D1 and D2 
and the temperature difference between the two points 
of the sphere in front of the detectors where the tem¬ 
perature T1 and T2 were measured (T1 and T2 were 
taken in front of D1 and D2, respectively). The count 
rate of the detector close to the higher temperature point 
was lower than the one of the detector close to the lower 
temperature point (the maximum temperature difference 
amounted to 0.05 K, while the maximum temperature 
excursion was from 13.5 to 17 °C). Finally, we were 
able to reproduce the count rate modulation by switch¬ 
ing on and off a heater placed at 1 m distance from the 
lead shield. 

Our conclusion is that large changes in the count 
rate in addition to the exponential decay are due to the 
displacement of radon and its progeny inside the glass 
sphere volume triggered by a small temperature gradient 
among different points on the sphere surface (as shown 
in Figure 2, a difference of 0.01 K is enough to produce 
a sizable effect on the count rate). The displacement of 
the radioactive atoms gives then rise to a different de¬ 
tection efficiency. 

3. Radon in liquid 

For a precise study of the ^^^Rn half-life it is nec¬ 
essary to suppress the temperature effect described in 
the previous section. Therefore, we selected an option 
able to ’immobilize’ radon and its progeny but allowing 
at the same time for a rather high radon concentration 
and, as a consequence, a high count rate. As a matter of 
fact, a statistical analysis shows that a precision on the 
life time better than one part over 10"^ requires an initial 
count rate of at least 400 cps and a measuring time of 
about 2 months. To achieve this we diffused radon in 
olive oil which has much higher viscosity than air and 
which permits in equilibrium conditions a radon con¬ 
centration 29 times higher than in air ifTbll . 

For the half-life measurement we minimized the 
background and its fluctuations by running the exper¬ 
iment underground at LNGS, where the muon and neu¬ 
tron flux are suppressed by six and three orders of mag¬ 
nitude, respectively, and by having a copper cube (10 
cm side, 2 mm thick walls) to contain the olive oil 
charged with radon. In particular, the copper cube filled 
with olive oil was charged with radon in the external 
laboratory by slowly pumping and diffusing through it 
the air coming from the cylinder containing the ura¬ 
nium rich rock. After about one week the cube was iso¬ 
lated, detached from the cylinder, hermetically closed 
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Figure 3: Spectra collected underground before charging the oil with 
radon (lower graph) and during a radon run (upper graph). 

and moved underground. The detector was a 3”x3” Nal 
powered and read out by a digiBASE (TM). The shield¬ 
ing to absorb the external gamma rays was made of at 
least 15 cm of lead and the laboratory temperature was 
kept between 12 and 13 °C. 
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Figure 4: Dead time as function of count rate for the half-life mea¬ 
surements. A detailed view at the beginning of the runs is given in the 
insert. Absolute deviations from a linear dependence are shown in the 
lower panel. 

Figure 3 shows the spectra collected underground be¬ 
fore charging the oil with radon and during a radon run. 
All the peaks due to the gamma rays given in Table 1 or 
to their sum are clearly visible. In addition, there is a 
small bump at 511 keV due to e^e “ annihilation from 
gamma rays interacting in the lead shield. In the analy¬ 
sis we consider the entire energy spectrum above 6 keV 
and not only the full energy peaks. This because we 
want to avoid any inaccuracy coming from the fitting 
procedure and we also want to increase the total rate in 


order to improve the statistics. 

Spectra are collected during consecutive time win¬ 
dows of 3600 s each, with the timing provided by the 
internal quartz oscillator of the digiBASE (TM). Its pre¬ 
cision and stability (better than 5 ppm/year) are enough 
for our purposes. The time stamp of each 3600 s mea¬ 
surement is then converted into the UTC (Coordinated 
Universal Time) time scale using a periodic check of the 
internal time of the data acquisition PC (free running) 
with respect to the UTC itself. This conversion gives 
rise to a relative systematic uncertainty on the life time 
measurement of less than 2-10^®. 

In the period from May 2014 to January 2015 we took 
four different measurements of the ^^^Rn activity to ex¬ 
tract the half-life and to measure its stability. Before 
charging with ^^^Rn we took a 20 day long background 
run with the copper box filled with olive oil and closed 
inside the lead shield underground. We took this mea¬ 
surement to preliminary check the background stability 
over time. 

Dead time, provided by the digiBASE (TM) with 20 
ms resolution, amounts to 0.2% at the beginning of the 
radon measurement, linearly decreasing with the count 
rate, as shown in Figure 4. From the lower panel in Fig. 
4 we see that, for a constant count rate, there is a maxi¬ 
mum relative fluctuation of the dead time of 4-10“^. As 
a consequence, the maximum relative systematic uncer¬ 
tainty on the life time measurement due to dead time 
fluctuations amounts to at most 8-10“®, one order of 
magnitude smaller than the statistical error. 



Figure 5: Count rate as function of time and normalized residuals for 
the second half-life measurement. The shaded area corresponds to the 
±2 <t uncertainty region. 

The count rate as function of time and the normalized 
residuals for the second measurement are shown in Fig¬ 
ure 5 (for the other three measurements the figures are 
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Table 2: The parameters from the fit A+B-e to the count rate as function of time. 


Run 

Running time [h] 

A [cps] 

B [cps] 

C[d] 

reduced 

1 

1301 

1.0209(9) 

452.98(5) 

5.51335(46) 

1.13 

2 

1462 

1.0239(7) 

457.98(5) 

5.51302(44) 

1.06 

3 

1185 

1.0247(10) 

463.16(5) 

5.51352(46) 

1.09 

4 

1357 

1.0241(9) 

354.82(4) 

5.51289(50) 

1.03 


essentially the same). For each measurement we ht the 
count rate, after dead time correction, with the function 
A+B e In Table 2 the values of the parameters and 
of the reduced chi squared are given, with computed 
by taking into account only statistical uncertainties from 
the Poisson distribution. 

The ht parameter A corresponds to the background 
which, as shown in Table 2, is slightly increasing dur¬ 
ing the measurements. Figure 6 shows the lowest en¬ 
ergy part of the background spectrum taken in one hour 
before starting the radon measurement and during the 
last hour of the third radon measurement, when radon 
has completely decayed. It is interesting to note that the 
46.5 keV peak due to ^'°Pb, which is the background 
component expected to increase with time during the 
radon measurement because of its long half life, has in¬ 
creased by less than 310“^ cps (the 46.5 keV peak in the 
background spectrum before hlling with radon is due to 
^'°Pb decay in the lead shield). If we impose into the ht 
a slope on the background of 10“^ cps/60 days we then 
have a relative variation of the half life of 7 10“^, one 
order of magnitude smaller than the statistical error. 



Figure 6: The lowest energy part of the background spectrum taken 
in one hour before starting the radon measurement (shaded area) and 
during the last hour of the third radon measurement (points). Sta¬ 
tistical uncertainties are shown only for the radon measurement. In 
addition to the 46.5 keV peak from ^*®Pb decay, the L (around 12 
keV) and K (around 75 keV) X-rays of lead are clearly visible. 


From the four measurements we obtain the follow¬ 
ing values for the half-life; 3.82157(32) d, 3.82134(30) 
d, 3.82169(32) d and 3.82124(35). They are perfectly 
compatible and, as a consequence, the ^^^Rn half-life 
has been constant, within 8 parts over 10^ (1 sigma er¬ 
ror for each measurement) during 9 months from May 
2014 to January 2015. We are now planning to continue 
these half-life measurements for at least another year. 

Finally, if we average the results of the four in¬ 
dependent and mutually consistent measurements we 
obtain the most precise value for ^^^Rn half-life: 
3.82146(16)s,ar(4)jj,5, d, where the value of the system¬ 
atic uncertainty is the sum in quadrature of all the pre¬ 
viously described systematic uncertainties which are 
given in Table 3. Our value has to be compared to the 
most precise recent measurement of 3.8224(18) d, ob¬ 
tained with a ^^^Rn aqueous solution in liquid scintilla¬ 
tor ifTTl . and to the world average of 3.8232(8) d ifT^ . 
We observe that the world average is a weighted average 
of seven measurements m, several of which are very 
old, with limited description of the systematics but with 
very small error, much smaller than in IfTTl . 


Table 3: The relative systematic uncertainties on the ^^^Rn half-life 
(in 10“’’ units). 


Time conversion 

Dead time fluctuations 

TTUpb 

Total 

2 

8 

7 

10.8 


4. Conclusions 

We have studied the decay of ^^^Rn at the Gran 
Sasso Laboratory by detecting the gamma rays from its 
progeny. The purpose has been the search for time mod¬ 
ulations of radioactivity by precisely measuring several 
times per year the half-life of a short lifetime (days) 
source instead of measuring over a long period of time 
the activity of a long lifetime (tens or hundreds of years) 
source. As source we have used radon diffused into 
olive oil: this way we removed the large fluctuations 
in the count rate triggered by a temperature gradient at 
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the surface of the volume containing air charged with 
radon (a difference of 0.01 K is enough to produce a siz¬ 
able effect on the count rate). Our results do not support 
the hypothesis of a possible influence of a component 
in the solar irradiation affecting the nuclear decay rate. 
In particular, we performed four half-life measurements 
in the period from May 2014 to January 2015 which 
provide a constant ^^^Rn half-life with a 8-10“^ preci¬ 
sion. Finally, we improve by one order of magnitude 
the precision on ^^^Rn half-life providing the value of 
3.82146(16),,d. 
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